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The frustrated pyrochlore antiferromagnet Gd 2 Ti 207 has an unusual partially-ordered magnetic 
structure at the lowest measurable temperatures. This structure is currently believed to involve four 
magnetic propagation vectors k£(|||)*ina cubic 4-k structure, based on analysis of magnetic 
diffuse-scattering data [J. Phys.: Condens. Matter 16, L321 (2004)]. Here, we present three pieces 
of evidence against the 4-k structure. First, we report single-crystal neutron-diffraction measure¬ 
ments as a function of applied magnetic field, which are consistent with the selective held-induced 
population of non-cubic magnetic domains. Second, we present evidence from high-resolution pow¬ 
der neutron-diffraction measurements that rhombohedral strains exist within magnetic domains, 
which may be generated by magneto-elastic coupling only for the alternative Tk structure. Finally, 
we show that the argument previously used to rule out the 1-k structure is flawed, and demonstrate 
that magnetic diffuse-scattering data can actually be htted quantitatively by a 1-k structure in 
which spin fluctuations on ordered and disordered magnetic sites are strongly coupled. Our results 
provide an experimental foundation on which to base theoretical descriptions of partially-ordered 
states. 

PACS numbers: 75.10.Jm,61.43.Bn,75.25.-j 


Geometrical frustration is a central theme of 
condensed-matter physics because it generates exotic 
magnetic states. These states can usually be divided into 
two categories: spin liquids, in which frustration inhibits 
long-range magnetic order, and spin solids, in which per¬ 
turbations to the dominant frustrated interactions even¬ 
tually generate magnetic order [T]. Defying this classifi¬ 
cation, some frustrated magnets exhibit partial magnetic 
order HHH], in which only some spins order at low tem¬ 
perature. Partial magnetic order is surprisingly preva¬ 
lent in frustrated systems [ng and partial structural 
order plays a key role in determining the properties of 
materials such as fast-ion conductors and high-pressure 
elemental phases inmn]. The motivation for studying 
partial order lies in understanding the interplay between 
order and disorder that may drive its formation [llHlJj . 
To achieve this aim, experimental determination of the 
nature of partial order in a real material is essential. 

The frustrated antiferromagnet Gd2Ti207 is a canon¬ 
ical partially-ordered system in which magnetic Gd^"*" 
ions (S' = 7/2) occupy a pyrochlore network of corner¬ 
sharing tetrahedra M- Two magnetic phase transitions 
occur at Ti = 1.1 K and T2 = 0.75 K to differ¬ 

ent states with magnetic propagation vector k = [555]* 
EHl HU- Here, we consider only the low-temperature 
(T ^ T2) state. Two different magnetic structures are 
equally consistent with most experimental data, includ¬ 
ing local probes [m HI ng and the magnetic Bragg 
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FIG. 1: (a) 1-k and (b) 4-k magnetic structure candidates 

for Gd 2 Ti 207 . In both structures, the ordered site is shown in 
red, the disordered site is shown in grey, and average spin ori¬ 
entations are shown as red arrows. A single crystallographic 
unit cell is shown; spin orientations are reversed in adjacent 
unit cells. 


scattering observed in powder neutron-diffraction exper¬ 
iments US HU- These candidate magnetic structures are 
shown in Fig. The “ 1 -k structure” involves a single 
k = [555]* while the “ 4 -k structure” involves a super¬ 
position of the four symmetry-equivalent ks (555)*- In 
both structures, 3/4 of the spins participate in long-range 
order, while the remaining spins show mostly short-range 
(paramagnetic) correlations with only a small ordered 
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magnetic moment m- Yet, the nature of partial order 
is very different for the two structures: the 1 -k struc¬ 
ture has rhombohedral symmetry (magnetic space group 
Ri3m) and disordered spins are separated by 7.2 A, while 
the 4-k structure has cubic symmetry (magnetic space 
group Fs43to) and disordered spins are separated by 
3.6 A. Evidence for the magnetic structure of Gd 2 Ti 207 
has so far been indirect: it was argued in Ref. [21] that the 
magnetic diffuse scattering observed in powder neutron- 
diffraction experiments is consistent only with the smaller 
distance between disordered spins in the 4-k structure. 

In this Letter, we present experimental evidence that 
the low-temperature magnetic structure of Gd 2 Ti 207 is 
1-k, challenging the 4-k structure that was previously 
proposed [5T]. We employ a combination of single-crystal 
and powder neutron diffraction to demonstrate the ex¬ 
istence of magnetic domains with non-cubic symmetry, 
consistent with the 1-k structure. We further show that 
ordered and disordered sites are not independent as was 
previously assumed HU, but are actually coupled via spin 
fluctuations. 

Direct evidence against the 4-k structure is obtained 
from single-crystal neutron-diffraction measurements. 
Magnetic domains of the 1-k structure select different 
k G {\W)* therefore have different single-crystal 
diffraction patterns, while 4-k domains are related only 
by translational and time-reversal symmetries and hence 
appear identical to neutrons. However, a macroscopic 1- 
k sample is expected to contain equal populations of the 
four degenerate k G ( 555 )* domains, which would yield 
identical magnetic Bragg scattering to the 4-k struc¬ 
ture. We addressed this problem by applying a weak 
magnetic field to our sample in order to break the do¬ 
main degeneracy. The WISH diffractometer [23] at the 
ISIS neutron source was used to measure magnetic Bragg 
peaks at T = 0.07 K with a magnetic field applied along 
the [lIO] direction. The sample was a single crystal of 
volume ^ 10 mm^, which was cut from a larger crystal 
of 99.4% isotopically-enriched ^®°Gd 2 Ti 207 prepared by 
the floating-zone image furnace method [251 [ 25 ] using an 
isotopically-natural Gd 2 Ti 207 seed crystal. 

The field-dependence of selected magnetic Bragg peaks 
is shown in Fig. Magnetic Bragg peaks in the {hhl) 
plane disappear on the application of a magnetic field, 
while the intensity of magnetic Bragg peaks outside the 
{hhl) plane increases. This result is consistent with the 
1 -k structure, because the applied field makes the same 
angle (90°) with two of the k which lie within the {hhl) 
plane, and a different angle (35°) with the remaining two 
k which lie outside the {hhl) plane (see SI). Hence one 
pair of domains is populated by the application of a field 
and the other pair is depopulated. The only scenario in 
which our data might be consistent with the 4-k struc¬ 
ture is if the applied magnetic field actually caused a 
magnetic phase transition rather than a domain imbal¬ 
ance. This is unlikely, because there is no experimen- 
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FIG. 2: Intensities of selected single-crystal magnetic Bragg 
peaks at different values of the applied magnetic field poH 
parallel to [110] at T = 0.07 K. Magnetic Bragg peaks are 
labelled in each panel; left-hand panels show peaks within 
the {hhl) plane and right-hand panels show peaks outside 
the {hhl) plane. Points are coloured as follows: /roH = 0 
(black hollow circles), poH = 0.1 T (blue filled diamonds), 
poH = 0.2T (green filled triangles), and poH = 0.5T (red 
hollow squares). For clarity, only /roH = 0 and poH = 0.5T 
are shown in the right-hand panels. 


tal evidence for a field-induced phase transition in either 
specific heat [18] or torque magnetometry [27] measure¬ 
ments for fields along ( 110 ) of less than 2 T, whereas we 
observe changes in peak intensity for much smaller fields 
B « 0.2 T. However, the evidence for field-induced spin 
reorientations in Er 2 Ti 207 |28j motivated further exper¬ 
imental work, which we discuss below. 

We now show that a signature of non-cubic domain 
structure persists in zero applied field. If the magnetic 
structure is 1 -k, magneto-elastic coupling may generate a 
rhombohedral lattice distortion (see, e.g., [ 25 ]), whereas 
no such distortion is possible for the 4-k structure. To 
look for a distortion, we measured the powder neutron- 
diffraction pattern of Gd 2 Ti 207 at T = 1.1 K and T = 
0.03 K. The HRPD diffractometer at ISIS was used be¬ 
cause it has among the highest-available reciprocal-space 
resolution {AQ/Q « 5 x 10“^). The powder sample was 
the same as used in previous studies [201 dH and was 
mounted in a Cu holder to which deuterated isopropyl 
alcohol (d-IPA) was added to improve thermal contact. 
A dilution refrigerator was used to reach sub-Kelvin tem¬ 
peratures. We performed a series of Rietveld refinements 
to the HPRD data in which the extent of rhombohedral 
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distortion parameter, D 
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FIG. 3: (a) Dependence of goodness-of-fit 7?wp on the dis¬ 

tortion parameter D (defined in the text) for Rietveld re¬ 
finement to high-resolution powder neutron-diffraction data. 
Green diamonds show refinement of the nuclear phase to 
T = 1.1K(> Ti) data, blue squares refinement of the mag- 
netic-l-nuclear phase to T — 0.03 K(<C T 2 ) data, and red 
circles refinement of the nuclear phase only to T = 0.03 K 
data, (b) High-resolution powder neutron-diffraction patterns 
at T = 0.03 K. Experimental data are shown as black circles, 
Rietveld fit as a red line, and data-fit as a blue line. The up¬ 
per line of tick marks indicates the positions of nuclear Bragg 
peaks and the lower line of tick marks indicates the positions 
of magnetic Bragg peaks. The region around 100 ms was ex¬ 
cluded as it contains a large background contribution from 
d-IPA that could not be robustly refined. 

distortion was systematically varied. This was achieved 
by defining a distortion parameter, D = Ch/'s/fioh — 1, 
where Ch and Ch are lattice parameters of the crystallo¬ 
graphic unit cell in the hexagonal setting of space group 
R3m. For the undistorted cubic structure, D — 0. Ri¬ 
etveld refinements were performed using Topas Aca¬ 
demic (version 5) [50] and the Isodistort software was 
used to generate structural models [5I1|51|. The Bragg 
peaks arising from Cu and V were fitted by Pawley re¬ 
finement and the unit-cell volume, sample peak-shape, 
and background parameters were refined. 

The dependence of the goodness-of-fit (described by 
the weighted-profile i?-factor i?wp) on D is shown in 
Fig. I^a). For refinement to the T = 1.1 K data, i?wp bas 
a symmetric minimum around H = 0 ; hence no lattice 
distortion exists in the paramagnetic phase within the 
resolution of our data. For refinement to the T = 0.03 K 
data, we fixed the magnetic structure to be 1 -k (we note 
that this does not bias the results, since the 1-k and 4- 


k structures have identical powder-averaged Bragg scat¬ 
tering for H = 0 [20]). We allowed two additional pa¬ 
rameters to refine to fit the magnetic phase: the size 
of the ordered magnetic moment /iord for the ordered 
3/4 of the spins, and a Gaussian broadening parame¬ 
ter for the magnetic peaks that corresponds to a finite 
magnetic-domain size ^domain- Our key result is that the 
minimum i?wp is now obtained for a non-zero rhombo- 
hedral distortion D = 0.00022(4) with lattice parameters 
Oh = 7.19352(9) A and Ch = 17.6244(6) A. The fit-to-data 
for this model is shown in Fig. ib). The rhombohedral 
distortion is statistically-significant: |ch — c[J|/(T(ch) = 4, 
where cl = 17.62177(3) A is obtained for refinement to 
the T = 0.03 K data with the cubic constraint that 
D — 0. We find /^ord = 6.7(1) for the ordered site, 
consistent with previous studies [iniiii]. Finally, we in¬ 
vestigate whether the sensitivity to this distortion arises 
primarily from the magnetic or nuclear peaks, by per¬ 
forming a further set of refinements to the T = 0.03 K 
data where the magnetic phase was not included in the 
model. We find that the sensitivity of the data to a rhom¬ 
bohedral distortion is greatly reduced for these refine¬ 
ments. This result suggests that the rhombohedral dis¬ 
tortion does not persist over the length-scale of the pow¬ 
der crystallites but only over the finite magnetic-domain 
size ^domain = 2.4(2) X 10^ A —an unusual behaviour that 
has also been reported for the small monoclinic distortion 
present in MnO at low temperature [53 [55] . 

Our observations of rhombohedral symmetry-breaking 
apparently contradict a previous analysis |21j . where 
the distances between disordered spins in the 1-k struc¬ 
ture were shown to be incompatible with the presence 
of a broad magnetic diffuse-scattering peak at Q « 
1.1 A“^. Implicit in this argument is the assumption 
that the ordered site does not give rise to diffuse scat¬ 
tering. While this assumption may appear plausible, 
it is actually flawed, as we now show. The integral 
of the magnetic scattering intensity over all Q and en¬ 
ergy, /int, is proportional to the total squared magnetic 
moment, -I- 1 ), which is temperature- 

independent [34] . For a site which shows long-range 
spin order, one component of spin is aligned along a lo¬ 
cal axis with maximum projection ±S at low tempera¬ 
ture, generating Bragg scattering with intensity propor¬ 
tional to The transverse spin com¬ 

ponents undergo fluctuations (e.g., spin waves) and con¬ 
tribute the rest of lint [SS] [36] . In a diffraction exper¬ 
iment, an energy integral over spin fluctuations is per¬ 
formed, so that magnetic diffuse scattering is present 
even in conventionally-ordered magnets at low temper¬ 
ature [571I5H]. In Gd 2 Ti 207 , we estimate the proportion 
of the diffuse scattering which arises from ordered and 
disordered sites by taking g = 7.94/tb for both sites 
and assuming that the ordered site has its maximum 
Mord = 7.0 /tb whereas the disordered site has /iord = 0. 
Then, 25% of Imt will be diffuse scattering from the dis- 
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ordered site, and 3/[4(S'+l)] = 17% of lint will be diffuse 
scattering from the ordered site. Hence, the net contribu¬ 
tion to the diffuse scattering from ordered and disordered 
sites is actually expected to be similar. 

We now develop a quantitative model of the partial- 
order state which is consistent with all our previous re¬ 
sults. We use reverse Monte Carlo (RMC) refinement 
[371 - 140] to determine whether spin configurations exist 
that obey three constraints, (i) A single model is consis¬ 
tent with the total (Bragg-|-diffuse) magnetic neutron¬ 
scattering data, (ii) Each Gd^’*' ion has the same fi. 
(iii) The component of the magnetic structure that shows 
long-range order (“average structure”) is 1-k. Spin con¬ 
figurations were refined to the powder-diffraction data of 
Ref. HU which were collected at T = 0.25 K using the 
D20 diffractometer at the Institut Laue-Langevin m- 
Nuclear Bragg peaks and non-magnetic background were 
removed from the data before refinement; full details of 
the data processing and RMC refinement algorithm are 
given in SI. 

The RMC fit is shown in Fig.|^a). The excellent qual¬ 
ity of the fit proves that the magnetic diffuse scattering is 
consistent with the 1-k structure. Refinements were also 
performed where the average structure was constrained 
to be 4-k, which yielded similarly good fits (not shown). 
The diffuse scattering shows two distinct components. 
First, there is a broad feature centred at Q ~ 1.1 
which arises from antiferromagnetic spin correlations at 
r « 3.6A [21]; second, sharper features are also observed 
at the base of magnetic Bragg peaks. To understand the 
origin of these features we calculate partial spin correla¬ 
tion functions, 


where p, cr label either ordered or disordered spins, Np is 
the number of spins of type p, Zpa.{r) is the number of 
spins of type a which coordinate a spin of type p at dis¬ 
tance r, and spins are normalised to unit length. We 
split the calculated (S(0) • S(r))^^ into two parts: the 
average-structure correlation function, which does not 
decay with distance and gives rise to Bragg scattering, 
and the spin-fluctuation correlation function, which de¬ 
cays with distance and gives rise to diffuse scattering. 
These quantities are shown in Fig. Bb). The magnitude 
of average-structure correlations is much greater for the 
ordered than the disordered site, as anticipated from the 
relative values of pord m- The spin fluctuations resem¬ 
ble the average structure and persist over several crys¬ 
tallographic unit cells, consistent with a superposition of 
low-energy spin waves |38|. 

While the intra-site spin correlations account for the 
relatively sharp diffuse scattering at the base of the mag¬ 
netic Bragg peaks, they are too long-ranged to account 
for the broad diffuse peak at Q « 1.1 A“^. Strikingly, 




FIG. 4: (a) Reverse Monte Carlo (RMC) fit to experimental 

magnetic total-scattering data from Ref. 1211 Experimental 
data are shown as black circles, the Bragg component of the 
RMC fit as a blue line, and the diffuse component of the 
RMC fit as a red line. A square-root scale is used on the 
vertical axis to highlight the weak diffuse scattering, (b) Par¬ 
tial spin correlation functions. Panels show correlations be¬ 
tween ordered-ordered (“0-0”), disordered-disordered (“D- 
D”), and ordered-disordered (“O-D”) sites. In each panel, 
average-structure correlations are shown as blue bars and 
spin-fluctuation correlations are shown as red bars. Note the 
different vertical scales in each panel. The cubic lattice pa¬ 
rameter a = 10.17 A. 


however, the inter-site correlation function shows a sig¬ 
nificant antiferromagnetic peak at r = 3.6 A. We there¬ 
fore propose that the broad feature arises from correlated 
spin fluctuations between the ordered and disordered sites 
of the 1 -k structure, rather than correlations within the 
disordered site of the 4-k structure as was previously 
proposed [21]. Importantly, this result shows that the 
ordered and disordered sites are not independent, as pre¬ 
vious studies suggested HHl HU , but are actually strongly 
coupled. However, this coupling can only occur via spin 
fluctuations, because symmetry constrains the inter-site 
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correlations of the average structure to be zero. While 
our neutron-diffraction data do not directly determine 
the timescale of these fluctuations, neutron spin echo 
[32] and muon-spin rotation [331 |33j measurements have 
shown that low-energy spin dynamics persist at T T 2 , 
suggesting that they are dynamic rather than frozen. 

Our determination of the nature of partial order in 
Gd 2 Ti 207 provides an experimental foundation for theo¬ 
retical studies which aim to understand the microscopic 
origins of partial-order states Our key result is 

that the magnetic structure of Gd 2 Ti 2 07 for T <C T 2 is 
not 4-k, but the 1-k structure is consistent with all our 
data. We note that our results apply only for T <C T 2 , 
and it remains possible that a different magnetic struc¬ 
ture could exist for T 2 < T < Ti, as was recently pro¬ 
posed in Ref. |T21 We anticipate that the coupling of 
disordered and ordered sites via spin fluctuations may 
prove important in directing the nature of the magnetic 
ground state. Inelastic neutron-scattering experiments 
could provide more detailed information on the nature 
of this coupling, but are extremely challenging due to 
the high neutron-absorption cross-section of isotopically- 
natural Gd. 
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